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^ ' Abstract 

CJ ' The Bekenstein-Hawking entropy suggests that thermodynamics is an intrinsic in- 

gredient of gravity. Here, we explore the idea that requirements of thermodynamic con- 
sistency could determine the gravitational entropy in other set-ups. We implement this 
idea in a simple model: static, spherically symmetric solutions to Einstein's equations 
corresponding to self-gravitating radiation. We find that the principle of maximum en- 
^ I tropy provides a consistent thermodynamic description of the system, only if the entropy 

OO ■ includes a contribution from the spacetime singularities that appear in the solutions of 

Einstein's equations. The form of the singularity entropy is stringently constrained from 
ff^ i consistency requirements, so that the existence of a simple expression satisfying these 

I constraints is highly non-trivial, and suggests of a fundamental origin. We find that the 

■ system is characterized by three equilibrium phases, and we conduct a preliminary in- 

vestigation of the associated phase transitions. These results demonstrate the point that 
gravitational entities other than horizons are endowed with thermodynamic properties. 

X 



^ ; 1 Introduction 

The current classical theory of gravity is very successful both in providing a broad theoretical 
framework and in its agreement with experiments. However, in its present form it is also 
incomplete in the sense that (i) it does not include in one integrated picture the causal notion 
of singularities; and (ii) it does not contain an intrinsic notion of entropy that it is of crucial 
importance for the study of black hole thermodynamics. The latter issue has motivated several 
ideas about the integration of gravity and thermodynamics. Some interesting examples are 
Jacobson's interpretation of Einstein's equations as equations of state [Ij, Padmanabhan's 
programme on gravity as a thermodynamic theory [2] and Verlinde's interpretation of gravity 
as an entropic force [3] — see also Ref. [1] for a related review. 

At present, most discussions on gravitational entropy focus on the entropy of horizons 
(black hole or particle horizons), for which precise mathematical expressions exist. However, 
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horizons are not the only gravitational entities to which entropy can be assigned. In particular, 
Penrose has proposed the association of entropy with gravitational singularities as a way to 
obtain a thermodynamically consistent description of cosmology [5]. A concrete expression for 
the entropy of singularities, which is currently missing (at least for singularities in Lorentzian 
space-times), would greatly increase the scope of the search for a general theory of gravitational 
entropy. 

The aim of this work is to present a concrete expression for the entropy associated to a 
class of spacetime singularities, using elementary principles of general relativity and thermo- 
dynamics. This result provides a direct verification of the conjecture that (at least some of) 
the singularities can be viewed as thermodynamic objects. 

The key idea is the following. If the gravitational degrees of freedom are carriers of en- 
tropy, then the total entropy in a gravitating system should be the sum of matter entropy and 
gravitational entropy (plus a small term due to their interaction.) However, the form of any en- 
tropy functional is constrained by requirements of thermodynamic consistency. In particular, 
a well-defined entropy functional should allow for the implementation of the maximum- entropy 
principle: in an equilibrium state, all unconstrained parameters describing a physical system 
take values that maximize the entropy functional. If one ignores the contribution of grav- 
itational entropy, the resulting entropy functional may fail to satisfy defining properties of 
entropy and thus, misidentify the physical equilibrium states. Hence, the requirement that 
the entropy functional satisfies basic thermodynamic consistency conditions could allow us 
to infer the form and properties of any contribution to entropy from the gravitational field. 
We show that this idea is successfully implemented in static, spherically symmetric solutions 
to Einstein equations describing self-gravitating radiation [HI [7] . It leads to an expression for 
gravitational entropy that is defined in terms of the properties the spacetime singularities that 
generically appear in these solutions. 

The key physical requirement in this paper is the consistency of the thermodynamic de- 
scription. It is, however, to be noted that the presence of a long-range force such as gravity 
leads to a significant modification of the usual thermodynamical properties — for reviews, see, 
Refs. [H [9l dni [11] . Most important is the loss of extensivity in the thermodynamic poten- 
tials. In ordinary equilibrium thermodynamics, the entropy is an extensive function on the 
system's thermodynamic state space F: S{XX) = XS{X), for any A > and X E T. This is 
not the case in gravity, even though in certain systems other scaling properties exist fT3] . 
The loss of extensivity implies that: the equilibrium configurations are inhomogeneous; the 
heat capacity may become negative, even though the specific heat of each volume element is 
positive [Ellis]; the microcanonical and the canonical distribution are not equivalent p^lTT] : 
and, crucially, matter entropy may not have global maxima for fixed values of energy |17] . 

Thus, many concepts of ordinary equilibrium thermodynamic concepts must be modified 
in order to fit the gravitational context. Nonetheless, black hole entropy suggests that the 
thermodynamic description constitutes a fundamental element of gravitational physics. For 
this reason, the fundamental principles of thermodynamics, as encapsulated in its basic laws, 
can be expected to hold even in a less stringent form. In this regard, we note that many con- 
cepts of the axiomatic approaches to equilibrium thermodynamics — for example [THl H^j — are 
applicable in the gravitational context. In these approaches, entropy is viewed as a function on 
a state space of macroscopic observables that determines the allowed transitions of a physical 
system, and thus it can be defined irrespective of whether the extensivity property is satisfied 
or not. In particular, the maximum-entropy principle, which is the key thermodynamic prin- 
ciple employed in this paper, is based on the statistical interpretation of entropy and on the 
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second law of thermodynamics. These are universal principles and they are expected to apply 
also to gravitating systems. 

With these considerations in mind we revisit the thermodynamics of static, spherically 
symmetric solutions to Einstein equations corresponding to radiation in a box. The station- 
arity assumption allows us to employ notions of equilibrium thermodynamics — entropy in the 
full gravity theory would require notions of non-equilibrium thermodynamics. Spherical sym- 
metry implies that the local degrees of freedom of the gravitational field are frozen, so that 
any contribution to the gravitational entropy would be 'topological'. In the present context, 
this means that gravitational entropy can be identified with the entropy of singularities — or 
possibly with entropy associated to internal boundaries. We chose the case of self-gravitating 
radiation [6], because its thermodynamic description is simple: the equation of state is linear 
and it contains no scale parameters. This allows for a semi-analytic treatment of the system 
and a relatively straightforward determination of the gravitational contributions to entropy. 
Furthermore, the equation of state for radiation is expected to be physically meaningful even 
when the energy density varies over many orders of magnitude. 

The basic thermodynamic variables for self-gravitating radiation in a box are the ones 
accessible to an observer outside the box, namely, the area AttR"^ of the box, the local temper- 
ature T of the box and the Arnowitt-Deser-Misner (ADM) mass M. The variables R, M and 
T define the boundary conditions of Einstein' equations for this system. They can be varied 
independently, each choice of {R, M, T) defining a different solution to Einstein's equations. 
However, R, M and T are not independent, when viewed as thermodynamic variables. For 
fixed size of the bounding box, the temperature of radiation is a function of the total energy. 
Thus, states of thermal equilibrium satisfy a structure relation of the form f{R,M,T) = 0, 
which cannot be determined by Einstein's equations alone. The standard treatment of static 
spherically symmetric solutions to Einstein's equations introduces an additional assumption, 
of regularity at the center (see, Ref. [6l [7] for self-gravitating radiation). Regular solutions 
form a set of measure zero in the space of all solutions; most solutions have a conical singular- 
ity at the center. The restriction to regular solutions follows mainly from the desire to avoid 
gravitational singularities, but it also provides the additional condition that is necessary in 
order to derive the thermodynamic structure equation f{R,M,T) = 0. 

However, the regularity condition has no thermodynamic justification. For example, it 
is not obtained from the extremization of any thermodynamic potential. If we assume that 
general relativity arises as the macroscopic limit of an underlying microscopic theory, the 
lack of a fundamental justification for the regularity condition poses a severe problem. There 
should be some physical mechanism guaranteeing the stability of regular solutions against 
statistical or quantum fluctuations of the microscopic degrees of freedom. At the macroscopic 
level, such a mechanism should be described in the language of thermodynamics. 

Regular solutions exist only in a specific sector of the thermodynamic state space. Outside 
this sector, the regularity condition cannot identify the equilibrium states. Again, this is 
problematic from the thermodynamic point of view. To see this, we note that for regular 
solutions there is an upper bound Tmax{R) on the temperature of the bounding box for a 
given value of i? [SI E] • Suppose that we bring the box into contact with a thermal reservoir 
at temperature T > T^axiR)- One expects that the box also acquires temperature equal to T 
|21j . But there exists no regular solution with this property, so we have no way to determine 
the resulting equilibrium configuration. Thus, the regularity condition does not suffice to 
describe even elementary thermodynamic operations to the system. 

In ordinary equilibrium thermodynamics, when the external constraints (boundary con- 



3 



ditions) to a system do not suffice to determine the system's state, the maximum-entropy 
principle is invoked |20]. However, the entropy Srad of self-gravitating radiation does not have 
global maxima for fixed values of energy. Thus, the maximum-entropy principle does not 
recover the regular solutions, even at the limit of weak gravity. This is a severe consistency 
problem, because at this limit we expect the ordinary thermodynamics of a photon gas to 
provide an excellent approximation. 

The problems above are resolved by the assumption that an entropy Ssing is associated to 
the singularities characterizing the solutions to Einstein's equations. Ssing is defined only by 
the asymptotic behavior of the fields near the singularity, and its form is determined by the 
requirement that the total entropy Stot = Srad + Ssing is maximized at the regular solutions. 
Thus, the regularity condition arises as a consequence of the maximum- entropy principle. The 
maximum-entropy principle also determines the equilibrium configurations in the sector where 
no regular solutions exist. We find two phases of singular equilibrium solutions in addition to 
the phase of regular solutions. 

It is important to emphasize that there are strong constraints to the form of Ssing, arising 
from requirements of thermodynamic and mathematical consistency. These constraints are 
so stringent, that it is a priori unlikely that they can be satisfied accidentally, or even by 
fine-tuning. Hence, the existence of an expression for Ssing that satisfies the constraints is 
highly remarkable. It is even more so, because Ssing has a very simple form, and there 
are plausible arguments suggesting that it is unique. For this reason, we believe that the 
expression for Ssing, determined here, is of a fundamental origin and it extracts information 
that is deeply embedded into the structure of Einstein's equations. Further elaboration of 
these ideas into other systems could provide significant novel information towards a general 
theory of gravitational entropy. 

The association of an entropy Ssing to the singularities does not necessarily imply a com- 
mitment into the physical existence of singularities. As we explain in more detail in Sec. 3.2, 
we view Ssing as a mathematical term, defined by the asymptotic field values on the singu- 
larity, that is included in the total entropy in order to obtain a consistent thermodynamic 
description of the system. The expressions "entropy of singularity" and "singularity entropy" 
employed in this paper are to be understood in this sense. They do not imply the existence 
of the singularity as a physical entity that possesses entropy. 

The definition and properties of the entropy Ssing is consistent with the idea of an intrinsic 
thermodynamic character of gravity. Borrowing Padmanabhan's analogy of spacetime to an 
elastic solid [2l|, the singularities can be viewed as analogous to defects in the solid. For weak 
external stress, the solid's equilibrium state is smooth and regular, but if the stress increases 
beyond a critical value, no globally regular configuration exists: configurations with defects 
are entropically favored. The thermodynamic picture of gravity provides an important shift 
of perspective: the singularities are not pathologies in the description of the system, but 
they may indicate a different thermodynamic 'phase' that is entropically favored for certain 
boundary conditions. 

The structure of this paper is the following. In Sec. 2, we study in detail the properties of 
the singular solutions of the Oppenheimer-Volkoff equations for self-gravitating radiation in a 
spherical box. In Sec. 3, we define the singularity entropy Ssing and we show that the regular 
solutions are indeed obtained by a maximum-entropy principle. We also identify the phase of 
singular equilibrium solutions and study the corresponding phase transitions. In Sec. 4 we 
discuss our results and their implications. 
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2 Solutions of the Oppenheimer-VolkofF equations for 
self-gravitating radiation 

2.1 Structure equations and boundary conditions 

The thermodynamic system under consideration is a spherical box of area A = ttR^ containing 
a self-gravitating photon gas in thermal equilibrium. The local temperature T of the photon 
gas is determined by the energy density p, 

p = (1) 

where a is the Stefan-Boltzmann constant. The equation of state is P = |p, where P is the 
pressure. 

In what follows, we work entirely within the confines of the classical theory, so h does not 
appear. It is therefore convenient to choose units where a = 1. 

Assuming that the system is in equilibrium, we consider only static solutions to Einstein's 
equations. We denote the associated timelike Killing field as ^. Outside the box, the spacetime 
metric is the Schwarzschild solution with ADM mass M. Inside the box, the metric is of the 
form 



in the usual spherical-symmetric coordinates {t,r,6,(j)). The mass function m(r) satisfies the 
equation 

dm 9 , , 

- = (3) 

with the boundary condition m{R) = M. The energy density p(r) is obtained by the solution 
of the Oppenheimer-Volkoff equation 



dp _ 4p{m + Ivrr^p) 

~ 1-212 



(4) 



An external observer has access only to variables that can be measured outside the box, 
namely, the box's area A, its local temperature T and the ADM mass M. (The mass is 
determined by measurements of the acceleration of free-falling test particles near the box). 
By Eq. ([1]), the temperature T determines the value of the density at the boundary p{R). 
Thus, the knowledge of M = m{R) and T = [p(i?)]^/^ allows us to integrate the system of 
Eqs. ([3]) and (jlj) from the boundary r = R towards the center r = 0, and thus identify the 
spacetime geometry inside the box. 

The system of Eqs. ([3H1]), integrated from the boundary r = R inwards, admits two 
different classes of solutions, distinguished by their behavior as r — )■ 0. The generic case 
corresponds to solutions with a conical singularity at r = 0. In these solutions, m(0) = — Mq 
for some positive constant Mq, and p ~ as r — i- 0. The other class consists of regular 
solutions; these satisfy m(0) = and p(0) = pc for some constant pc > 0. The regular 
solutions form a set of measure zero in the space of all solutions. They define spacetimes that 
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are everywhere locally Minkowskian. In Sec. 2.3, we show that no solutions with m(0) > 
exist when Eqs. ([3]) and (jlj) are integrated from the boundary inwards. 

In regular solutions, m(r) is positive valued and m(r — >■ oo) = M, hence, it is usually 
interpreted as the total energy within a spherical surface of radius r. This interpretation is 
heuristic, because in general relativity there is no invariant way of defining the total energy in a 
finite region. In singular, solutions, the function m(r) takes negative values near the singularity 
so its interpretation as total energy is untenable. In particular, the value m(0) = —Mq is a 
topological feature of the singularity, with no a priori relation to energy. As we explain in 
Sec. 3.3, a better measure of energy in singular solutions is obtained in terms of the Komar 
integral. 

Energy and temperature are not independent variables. In an ordinary (i.e., non-gravitating) 
photon gas the entropy 5* is a function of the energy U and the volume V, and the temperature 
T is defined as T^^ = {dS/dU)v ■ Similarly, in self-gravitating radiation, there are only two 
independent thermodynamic variables (for example, M and R). However, a unique solution 
of Eqs. ([3]) and (jl]) requires the specification of three independent variables, M, T and R. 

In usual studies of spherical symmetric solutions to Einstein's equations, one restricts at- 
tention to the subspace of regular solutions. For self-gravitating radiation, this implies solving 
Eqs. (|3]) and as a boundary- value problem, with m(0) =0 and m{R) = M. Regularity is 
an additional condition that determines the functional dependence of the temperature T on 
the variables R and M. 

The equation of state for radiation does not contain any length scale. For this reason, 
the Oppenheimer-Volkoff equation possesses a scaling symmetry: the transformation M — )■ 
Xm,R^ XR,T ^ \-^/^T preserves the solution curves of the system. This symmetry implies 
that the structure equations are equivalent to a two-dimensional time-homogeneous dynamical 
system. The definition of the variables 

e = In^ (5) 
2m(r) 

u = (6) 

r 

V = 47rr^p (7) 

brings Eqs. and into the form 

u' = 2v-u (8) 

v' = ^, 9) 

1 — u 

where the prime refers to differentiation with respect to C,- Eqs. ([8]) and (Q are integrated 
from ^ = (boundary, r = i?) to ^ — )■ — oo (center, r = 0). In what follows, we denote the 
values of u and v at the boundary as uq := u(^ = 0) and vq := v{C, = 0). Thus, uq = 2M/R 
and vo = 4:7rR^T\ 

2.2 Regular solutions 

Regular solutions to Eqs. ([8]) and Qj have been studied in detail in Refs. O [7] . In a regular 
solution, the functions u and v vanish at the center. This means that the set of regular 
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Figure 1: The curve on the uq — vq plane representing regular solutions to Eqs. dHHS]). 
solutions corresponds to the solution curve of the differential equation 



with initial condition v{0) = 0. This curve is plotted in Fig. 1. Each point in this curve 
corresponds to boundary data {uo,vo) that define a regular solution. 

Fig. 1 shows that regular solutions are not defined for all values of Uq or vq. The point 
P corresponds to the maximum value of Vq allowed by a regular solution, and the point Q to 
the maximum value of Uq allowed by a regular solution. We denote the coordinates of P as 
{up,vp) and those of Q as {uq,vq). Thus, for regular solutions Uq < uq and vq < vp. We 
find numerically that up ~ 0.3861, vp ~ 0.3416, uq ^ 0.4926 and vq ~ 0.2463. In terms of 
the variables {R, M, T), the inequalities above imply that for regular solutions, TR^/"^ < 0.406 
and MjR < 0.246. Hence, the restriction to the subset of regular solutions implies an upper 
bound to both temperature T and ADM mass M for a given area AnR"^ of the bounding box. 

Since the curve in Fig. 1 develops a spiral for uq > up, there are several different values 
of Vq corresponding to each value of uq, when uq G {up,uq). This is because the point P 
in Fig. 1 signifies the onset of thermodynamic instability: {dM/dT)p becomes negative at 
P [22|. However, restricting to the segment from the origin to the point Q, we can define a 
single-valued map Vreg{uo) that assigns a unique value Vq = Vreg{uo) to each uq G [0,uq]. 

2.3 Singular solutions 
2.3.1 A 'no-horizon' theorem 

In static spherically symmetric solutions to Einstein's equations an event horizon is identified 
by the condition 2m{r)/r = 1, or equivalently m(^) = 1. The singular solutions to Eqs. ([H]) and 
([9]) exhibit no event horizon, i.e., they are 'naked'. The integration of the Oppenheimer-Volkoff 
equation from the boundary inwards encounters no horizons. This statement is a special case 
of a general theorem valid for a larger class of equations of state p3], so the lack of horizons 
is not a special property of self-gravitating radiation. 

The theorem of Ref. [22] applies to solutions with continuous first derivatives of the energy 
density p(r). However, it is conceivable that the energy density p could be discontinuous or 



dv _ 2v{l - 2m - Iv) 



(10) 



du (1 — u){2v — u) 
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divergent on a horizon, a possibility that is not excluded by the theorem above. For this 
reason, we provide an elementary proof that no horizon is encountered in the integration of 
Eqs. ([3H1]), in which the assumptions about the behavior of p(r) on the horizon are weaker 
than those of Ref. [23]. The proof proceeds as follows. 

Assume that an horizon forms for a value ^ = ^* G (— oo,0), i.e., = 1. The horizon 

is reached by integration from the boundary inwards, and we are interested in the behavior 
of the solution as ^ — > For ^ > -u < 1, so we define e(^) := 1 — m(^) > 0. We assume 
continuity of first derivatives of u and v for all ^ G (— oo,^*), but make no such assumption 
for the behavior of u{C,) and v{C,) at ,^ = C,*+ There are three possibilities for the behavior of 
f (^) at the limit ^ — t- 

1. lim^_^^^^v{^) = 0. Then, in a sufficiently small neighborhood of .^o; Eq. ffTOj) is approxi- 
mated by 



dv 2v 
de e 



which is solved by f = ce~^ for a constant c. On the horizon (e = 0) f (^) — )■ oo, thus 
contradicting the initial assumption. The only exception is the trivial case c = 0, which 
corresponds to f (^) = for all ^ > 

lim^^^^^v{^) = V > 0. In a sufficiently small neighborhood of ^o, and assuming v ^ ^, 
Eq. (fTOj) is approximated by 



dz 2v{l + p 



de e{2v - 1) 



(12) 



where we wrote z = v—v. The solution is 2; = In ce 2^-1 for a constant c > 0. Then 



2i;(l+|_ii) 



on the horizon either 2;— j-ooor^;— 00, thus contradicting the initial assumption. 
Similarly, one can show that the case ^ = | also leads to contradiction. 

3. lim^^^^^v{^) = (X). In a sufficiently small neighborhood of Eq. flTUl) is approximated 
by 

dv 2v , , 

with solution v = ce^/^. At the horizon, w — )• 0, thus contradicting the initial assumption. 



We conclude that no singularity forms from the integration of Eqs. (I8H9]) from the boundary 
inwards, except for the trivial case that v{^) = for all > This exception corresponds 
to a boundary condition p{R) = 0, or equivalently T{R) = 0. Given an non-zero value of the 
ADM mass M, this solution represents a box containing a Schwarzschild horizon at r5 = 2M 
and no radiation. The solution cannot be continued to r < r^, because the Killing field ^ is 
spacelike in the Schwarzschild interior. 

A consequence of the no- horizon theorem above is that no solutions with m(0) > exist. 
For, if such a solution exists, it satisfies lim^^ow = 00 as r — )■ 0. Since -u < 1 at the boundary, 
continuity of u implies that u{r^) = 1 for some r^, G (0,i?), thus contradicting the no-horizon 
theorem. 
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The lack of an horizon in the singular solutions is essential for the consistent thermody- 
namic description of the system. This can be seen from the following argument. Let the 
system be initially in a state corresponding to a regular solution. We use a quasi-static, adia- 
batic process to compress the gas and bring it into a state corresponding to a value of Uq for 
which no regular solution exists. We then reverse this process and return the system to its 
initial state. If the singular solution involves a horizon, the reversed process would lead to the 
horizon's disappearance, and thus to a violation of the second law of black hole mechanics. It 
is therefore significant that the only solutions with an horizon have vanishing temperature at 
the boundary. These states cannot be reached in a finite number of steps by an initial state 
with T ^0. 

2.3.2 Properties of the singular solutions 

Numerical integration of Eqs. ([H]) and ([9]) from the boundary inwards shows that there exist 
two qualitatively different types of singular solutions, which we denote as type I and type II. 
Their main features are the following. 

Type I solutions correspond to small values of v near the boundary, and u increasing as 
we integrate from the boundary inwards. In the neighborhood of a point, Tc, u blows up, 
approaching, but never reaching, the value u = 1 that defines an horizon. Integrating the 
solution to r < rc, we find that u decreases rapidly and tends to — oo as r — )■ 0. The energy 
density p has a very sharp peak at a value of r slightly smaller than Tc, and then decays 
to at r = 0. This implies that most of the photon gas is concentrated in a narrow shell 
around r^. Since in this region u is very close to unity, these solutions describe a 'wall' formed 
around a 'deformed' horizon, similar to the brick wall configurations proposed by 't Hooft 
|24j . In type II solutions, u is an increasing function of r, and it tends to — oo at r — )■ 0. In 
these solutions, the density p typically exhibits a maximum at intermediate values of r. Some 
properties of type II solutions are discussed in Ref. [25]. 

The two classes of solutions above are not exhaustive, and they are defined in terms 
of qualitative characteristics rather than precise mathematical conditions. Nonetheless, the 
distinction is significant. For uq < up, type I solutions correspond to pairs (uo,fo) such that 
vq < Vreg{uo) and type II solutions to pairs {uq, vq) such that uq > Vreg{uo)- For up < uq < uq), 
the interchange between the two types of solution is rather complex, as there exist several 
regular solutions for each value of uq. For uo > uq the qualitative distinctions between the 
two types of solution are gradually lost. In general, for Vq/uq « 1 the solutions are of type 
I and for Vq/uq » 1 they are of type II. 

Fig. 2 plots the typical behavior of the mass function m(r) and the energy density p(r) as 
functions of r for the two types of singular solutions. The dimensionless variables m(r)/M, 
E?p{r)/M and r/R are employed in these plots. 

The two types of solution do not differ in their approach to the singularity at r = 0: 
m(0) = —Mo for some constant Mq > 0. From Eqs. (E]— H]), we find that for r near 

m(r) = -Mo + -kr^ (14) 
5 

p(r) = Ar^ (15) 

for some constant k. It is convenient to employ the dimensionless parameters po = Mq/ R and 
K = fci?^, because these are functions of uq and fo alone. This can be seen from Eqs. (IHHQ]), 
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Figure 2: The normalized mass function m{r)/M as a function of r/i? and the normahzed density 
E?p{r)/M as a function of r/i? for a type I solution (solid curve) and a type II solution (dashed 
curve). The sudden drop in m{r) at r ~ 0.05i? in the type I solution coincides with the peak in 
the density. Note that the peak in the density is so sharp that we had to use a logarithmic scale in 
the vertical axis, in order to include it in the plot. The type I solution plotted here corresponds to 
{uq,vq) = (0.25,0.22) and the type II solution to {uq,vo) = (0.25,0.40). 



whose solutions have an asymptotic behavior 

u = — 2/ioe~^ V = ne^^ (16) 

as ,^ — 7- — oo. 

The values of Mq and k are the same for all points (i?, M, T) of the same solution curve. 
This fact, together with the asymptotic behavior Eq. (fT6|) . imply that /^o and k satisfy the 
following equations. 

{2vo-uo}- \ ^ -— = -yUo (17) 

ouq 1 — Wo (yvo 

{2vo-uo)^ + ^ = 4:K. (18) 

The functions /^(mo, Vq) and k,{uo, Vq) will be important for the construction of the singular- 
ity entropy in Sec. 3. Numerical evaluation shows that they possess the following properties. 

1. For fixed Uq, Hq{uq^Vq) ~ Vq^ and k{uq,Vq) ~ Vq^ as Vq — )■ 0. 

2. For fixed Uq, /io('Uo,fo) ~ and k{uq,Vq) ~ Vq as Vq oo. 

3. For fixed mq, both /iq and diverge to infinity as vq approaches from below the value 
Vregiuo) characterizing a regular solution, n diverges more strongly: /Iq/k — )• as I'd — > 
Vregiuo)^. Both fuuctious are continuous in their approach to Vregiuo) from above. 
In particular, the function iXq/h is continuous at Vreg{uo) and is equal to zero. The 
characteristic behavior of /xq and k in the vicinity of Vregiuo) is illustrated in Fig. 3. 

In the vicinity of r = the metric Eq. ([2]) becomes 

ds^ = -(1 - uo)xl^-de + ^ + r\de^ + sin^ Odcj)''). (19) 
\ K r 2/io-R 
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Figure 3: Plots of In fj,o{uo,vo) and In k{uq,vq) for uq = 0.1 and different values of vq near the point 
^^0 = Vreg{uo) ~ 0.1379 Corresponding to the regular solution. 



The proper radius coordinate corresponding to the metric Eq. ( !T9|) is x = |r^/^/ y/2^QR. A 
two-sphere of proper radius x around r = has area equal to 4'7r(|/io-R)^'^^x'^''^. This impUes 
that the spacetime is not locally Minkowskian around r = (since in that case the area 
should be Attx'^), and that the point r = is a conical singularity. The fact that p — )■ as 
r — )■ implies that the stress-energy tensor vanishes on the singularity, and, hence, unlike 
conical singularities in two spatial dimensions |26], the conical singularities appearing in self- 
gravitating radiation cannot be interpreted as particles. 

A singularity not hidden within an event horizon is usually thought of as violating the 
cosmic censorship hypothesis (CCH) [27], and thus it may be deemed unphysical. The main 
motivation for the CCH is to guarantee the preservation of causality and predictability in 
general relativity, which are threatened by the existence of naked singularities. The greatest 
problem is posed by curvature singularities, such as the ones appearing in the black hole 
interior. However, the conical singularities appearing in self-gravitating radiation are rather 
'benign', in the sense that they do not give rise to inextensible causal geodesies. Causal radial 
geodesies satisfy the equation 

1 .2 /io-R , e^/io-R , - 
-r = -g + , 20 

2 r a 

where a = (1 — Uo)^^^/^R, e = ai/r is the energy corresponding to the geodesic, (7 = 1 for 
timelike geodesies and (yf = for null geodesies. The derivative in Eq. (!20|) corresponds to 
proper time in timelike geodesies and to an affine parameter A in null geodesies. 

From Eq. fl2Up it is evident that no timelike geodesies reach the singularity: incoming 

geodesies bounce back at a minimal radius Tmin = a/e^- Null geodesies satisfy f = ±t / ^^-H^ 



and they reach r = within a finite value of their affine parameter (and also of the coordinate 
time t). An incoming null geodesic reaching the singularity simply becomes an outcoming null 
geodesic. 
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3 Thermodynamic consistency and the entropy of sin- 
gularities 

3.1 The identification of singularity entropy 

As discussed in Sec. 2.1, self-gravitating radiation in a spherical box is described by two 
independent thermodynamic variables, which can be chosen to be the ADM mass M and the 
'radius' R. However, the structure equations ([3]) and (jlj) require the specification of three 
independent variables. 

In usual thermodynamic systems the values of any unconstrained parameters are deter- 
mined by the maximum-entropy principle, namely, the statement that their equilibrium values 
maximize the entropy for a given value of the energy [20]. In order to implement this princi- 
ple, it is necessary to construct an entropy functional that depends on the three independent 
variables required for the specification of a unique solution of the structure equations (|8]) and 
([9]). Thus, we must construct an entropy functional S{M,T,R) or equivalently S{uo,Vo, R). 

To this end, we first consider the total entropy Srad of the radiation in the box. The 
entropy density s of radiation is s = Hence, 



Srad = ^ dr ^ = T.(47r)^/M -==. (21) 

'0 /l _ 2rn 6 Jo y/l-U 



For solutions to Eqs. dHHO]), the integrand in Eq. fl2T]) is a total derivative [6], i.e., 

^1/2^3/4 ^ ^ ^ + 1^ X 

2 ^3/2 _ ^22) 



y/T^ dr \6vy^^/T^ 
Hence, the radiation entropy is a sum of two terms 

Srad{Uo,Vo,R) = Si{Uo,Vo,R) + So{Uo,Vo,R), (23) 

where 5*1 is defined by the field values at the boundary r = R 

S,iuo,vo,R) = ^(4vr)V^ ^^; + 1^ ^3/2 (24) 
and 5*0 by the field values at r = 

S„(u„, „„. R) = -?(4.)'/* li,„ (^;iii^r'/=) = 1(4.)'/* (itp^) fi3/._ 

9 r^O \vi/'iy/l - u / 3 V fi{Uo,Vo) J 

where the limit was taken using Eq. ( IT6|) for the behavior of the functions u and w as r — )■ 
— > — oo). We note that the term 5*0 vanishes for regular solutions. 

The usual study of static spherically symmetric solutions of Einstein's equations invokes 
the condition that the spacetime metric is regular at the center. The regularity condition is 
motivated by the perception of singularities as pathologies in the description of space-time 
geometry that should be avoided in order to guarantee the causality and predictability of the 
theory. However, as shown in Sec. 2.3.2, this is not a problem for the singular solutions 
considered in this paper. In Sec. 3.2, we will also argue that the consideration of singular 
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Figure 4: Srad/R^^"^ as a function of vq for different values of uq. The radiation entropy does not 
have a maximal value, it tends to infinity at the limits u o — )• and vq — )• oo. The dots denote the 
regular solutions. 

solutions does not necessitate the acceptance of singularities as genuine spacetime points, at 
least in the context of the approach developed in this paper. However, the a priori exclusion 
of singular solutions as unphysical, implies that there can be no description of the system 
in the sector of the state space where regular solutions do not exist. There will be total 
unpredictability for any thermodynamic operation that takes the system into this sector. 

Moreover, the lack of a more fundamental justification for the regularity condition is con- 
ceptually problematic. If general relativity arises as the macroscopic limit of an underlying 
microscopic theory, there should be some physical mechanism guaranteeing the stability of 
regular solutions against statistical or quantum fluctuations of the underlying structure. At 
the macroscopic level the only language we possess to describe such mechanisms is that of 
thermodynamics. 

There is therefore good reason for providing a thermodynamic interpretation of the regu- 
larity condition. This implies that one should make no a priori distinction between regular 
and singular solutions. The regular solutions should be identified as the physical states of the 
system by the fundamental condition of thermodynamic equilibrium, namely, the principle of 
maximum entropy. This means that regular solutions should correspond to global maxima of 
the entropy functional, for any given values of energy M and box area AttR'^. 

However, the radiation entropy Srad of Eq- ( 123]) does not allow for the implementation of 
the maximum-entropy principle. Srad has no global maximum for any fixed values of R and M. 
The typical dependence of Srad on vq for fixed values of Uq = 2M/ R and R is depicted in Fig. 
4. There, we see that Srad — >■ oo as f o — )■ and as — )■ oo. In fact, the regular solution is close 
to the global minimum of the entropy functional. This is quite problematic, because Srad has 
the same behavior even at the limit of near- vanishing gravitational interaction (i? >> 1 and 
M/R arbitrarily close to zero). This would imply that the ordinary phase of thermal radiation 
is thermodynamically unstable even in presence of a tiny gravitational self-interaction. 

The behavior of Srad, described above, is not unexpected: it is an analogue of the well 
known gravothermal catastrophe [17]. However, it means that the thermodynamic justification 
of the regularity condition cannot be obtained by the consideration of the radiation entropy 
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alone. The only solution to this problem is to assume that the entropy functional also includes 
a contribution from the gravitational field. In static spherical symmetric configurations, there 
are no local gravitational degrees of freedom, so one expects that any gravitational contribution 
to entropy is topological, i.e., it corresponds to the properties of the singularities (or, possibly, 
to the properties of an interior boundary). If the entropy of the singularities is significantly 
lower than the entropy of locally Minkowskian space, then the regularity condition could follow 
as a consequence of the maximization of entropy. 

A different line of thought suggesting that singularities in general relativity could be associ- 
ated to entropy originates from Penrose Penrose's idea is that different types of singularity 
are distinguished by their entropic content. In particular, strong curvature singularities such 
as the one in the Schwarzschild-Kruskal spacetime correspond to high entropy, while 'mild' sin- 
gularities (for example, the Big-Bang where the spacetime is conformally regular) correspond 
to low entropy. The primary aim of this proposal is to provide a consistent thermodynamic 
description of cosmology. However, the idea applies naturally to the present system. Loosely 
speaking, the singularities of self-gravitating radiation are mild (no geodesic incompleteness), 
so Penrose's conjecture would suggest that they correspond to low entropy, as would be re- 
quired for a thermodynamic explanation of the regularity conditions. 

Thus, the requirement of a consistent thermodynamic interpretation for the regularity 
condition leads to the conjecture that the total entropy of the system can is a sum of the 
radiation entropy Srad and a term Sgmg for the entropy of singularities. 

Stot = Srad + Ssing- (26) 

Regular solutions have maximum entropy, only if the entropy of the singularities is lower 
than the entropy of locally Minkowskian geometry. Otherwise, a locally Minkowskian space- 
time would be unstable even in absence of matter. It is convenient to choose Stot so that the 
solutions with M = have zero entropy for all R. The term Ssing is then negative-valued. 
It turns out that with this convention, the total entropy Stot may take negative values. This 
contradicts the physical requirement that entropy is positive-valued. To make Stot positive- 
valued, we must add a reference term Sfiat{R), corresponding to the gravitational entropy of 
fiat spacetime within a spherical box of radius R. However, Sfiat{R) does not depend on uq 
and vq, and hence, it does not affect the specification of the equilibrium states through the 
maximum-entropy principle. For this reason, Sfiat{R) is ignored in the present context. The 
entropy of fiat spacetime, viewed as a reference state for the system, cannot be fixed by purely 
thermodynamic arguments; it requires the knowledge of the microstates in the underlying 
theory. In Sec. 3.3, we show that the existence of a lower bound in Stot also implies a lower 
bound for Sfiat{R)- 

We point out that a non-zero value of Sfiat{R) does not contradict the idea that Minkowski 
spacetime is the unique lowest energy state of gravity, and thus, that it should be assigned 
zero entropy. Sfiat{R) refers to the entropy of the degrees of freedom localized within the 
box of radius R. Even if the entropy of a globally Minkowski spacetime is zero, a microscopic 
definition of Sfiat{R) would require tracing out the degrees of freedom outside the box. Hence, 
even if globally Minkowski spacetime is the vacuum in a quantum theory of gravity, the reduced 
state within the box is expected to be mixed due to correlations across the boundary, and thus 
possessing non-zero entropy. An analogous calculation for a scalar field leads to an entropy 
proportional to the area of the boundary [28] . 

Lacking a general theory for gravitational entropy, we can construct Ssing only by con- 
straining it with requirements of physical and mathematical consistency. These requirements 
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are the following: 

1. The application of the maximum-entropy principle to general self-gravitating systems 
ought to recover the successful theory of stellar stability, which is based on the regularity 
condition. This means that the regular solutions must correspond to global maxima of 
the entropy Stot- For self-gravitating radiation in particular, we should recover a large 
segment of the uq — vq curve of Fig. 1; at the very least the segment starting from the 
origin and ending at the point P. 

2. The radiation entropy S'rad has the following behavior under rescalings: Sradiuo,vo, ^R) = 
^^^"^Sradiuo, vq, R), for any A > 0. If Ssing has a different behavior under rescalings, then 
the value of Vq for maximum-entropy configurations would carry a dependence on R 
in addition to the dependence on uq, and we would not be able to recover the regular 
solutions as global maxima of entropy. This means that the singularity entropy should 
also rescale as Ssing{uo,vo, XR) = X^/'^Ssing{uo,vo, R). 

3. The existence of global maxima of entropy requires that Stot{uo,vo, R) is bounded from 
above, for fixed values of Uq and R. This implies, in particular, that Stot converges to 
finite values at the limits f o — )■ and — ?■ oo, and that neither of these values is a 
global maximum. For constant R, Stot should also be bounded from below, so that Stot 
can be made positive- valued by the addition of a term Sfiat{R) that has no dependence 
on uq and Vq. 

4. The term Ssing depends only on the properties of the singularity, i.e., on (limits of) 
values of the field variables at r = 0. Assuming that Ssing depends only on the values 
of m(r) and of its first derivative at the limit r — )■ 0, Ssing is a functional of Mq and k 
alone. 

The four conditions above are very restrictive and the existence of a functional Ssing that 
satisfies them is highly non-trivial. To find Ssing, we observe that Srad ~ Vq as f o — > and 
Srad ~ v^^'^ as fo — )■ oo. Ssing should havc the same asymptotic behavior with Srad if property 
3 is to be satisfied. The only combination of /iq and k, with this asymptotic behavior at the 
limit r — )■ 0, which vanishes in the absence of singularities, is (hI/kY^'^. 

Thus, given Eq. (l25l) . the natural candidate for the singularity entropy is a term Ssing = bSo 
for some constant b. We determine b by the asymptotic behavior of Srad/ Sq at the limits f o — > 
and f — !■ oo. Since two conditions are employed in order to determine a single parameter, the 
generic expectation is that the problem admits no solution. It is then quite remarkable that the 
two limits identify the same value of b. We find that limt,(,_5.o Srad/So = lim^g^oo Srad/ Sq = 6. 
Thus, 6 = — 6 is the only value compatible with the required asymptotic behavior of Ssing- It 
follows that 

Stot = Srad ~ 65*0 = 5*1 — 55*0. (27) 

The entropy Stot Eq- (l27|) satisfies conditions 2-4 above. The issue is whether Stot also 
satisfies condition 1, i.e., whether its global maxima correspond to regular solutions. Impres- 
sively, it does so. Stot/R^^'^ as a function of fo for different values of Uq is plotted in Fig. 5. 
For each uq, there is a clear maximum at a specific value of wq, which we denote as Veq{uo)- 
For uo < Uq, the maxima of Stot coincide with regular solutions. For vq < Veq{uo) the singular 
solutions are of type I, while for vq > Veqiuo) they are of type II. The dependence of Stot on 
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fo is different in the two types of solution. Tfiis results to non-smooth behavior of Stot at the 
maxima: the first derivative of Stot is discontinuous at vq = Veq{uo). 

The behavior of Stot changes for uq > up, as can be inferred from the presence of the spiral 
in Fig. 1. Stot develops additional local maxima, reflecting the presence of multiple regular 
solutions with same value of mq. These solutions are not maxima with respect to the radiation 
entropy Srad- The derivative of Stot is discontinuous at the local entropy maxima. The local 
maximum with the largest value of vo is also the global maximum, and corresponds to the 
equilibrium configuration. The local maxima correspond to metastable configurations. 

In general, metastable configurations in self-gravitating systems may be characterized by 
very long decay times |10j, and thus considered as physical. An estimation of the relevant 
timescales is not possible in the present context, because our analysis lies entirely within equi- 
librium thermodynamics; hence, we cannot make any statement about the physical relevance 
of the metastable solutions. However, we note that as uq increases, the local maxima ap- 
proach each other and at uq = uq they merge into a single maximum and the derivative of 
Stot becomes everywhere continuous. The existence of a single maximum persists for all values 
Uq> Uq. The fact that the distance between the local maxima becomes arbitrarily small near 
the critical point u = uq strongly suggests that the local maxima would strongly contribute 
in the properties of the phase transition in a statistical-mechanics description of this system. 

We also note that there is a large number of regular solutions for values of Uq near the 
center of the spiral in Fig. 1. We have verified that these regular solutions indeed correspond 
to local entropy maxima, at least within the degree of accuracy allowed by our numerical 
evaluation. 

Fig. 5 shows Stot for a small range of values of vq around the local maxima. However, an 
equilibrium solution corresponds to a global maximum of Stot- In Fig- 6, we plot Stot for a 
larger range of values of fo, capturing also its asymptotic behavior for small and large values 
of fo, and showing that the local entropy maxima are also global. Fig. 6 also illustrates the 
fact that the entropy Stot is bounded from below. 

The question arises whether the choice Sgmg = —QSq for the entropy of the singularities 
is the unique expression satisfying conditions 1-4 above. In general, if the four conditions 
admit one solution Ssmg, then any function f{Ssing) such that (i) f{Ssing) — ^ Ssmg at the 
limits fo — )■ and fo — )■ oo, and (ii) has the same global maxima with Sging, provides an 
alternative solution. In principle, there may exist many such functions, even though their 
explicit construction would be rather artificial. 

However, a physically natural assumption is that Sging is a boundary term obtained from 
the integration of a total divergence drS, where S a functional of the fields m(r) and p(r). It is 
also natural to expect that Sging depends only on the topological properties of the singularity 
which, given the scaling symmetry of the system, are the same in all geometries that correspond 
to the same solution curve of Eqs. ([H]) and This property is also satisfied by /io and k. 
Then, Ssmg should satisfy an equation of the same form as Eqs. (fT71 - [T8|) . namely, 

{2Vo - Uq)^-^ + ^ — ^ = -rSsing, (28) 

where 7 is a constant defined by the asymptotic behavior of 5 as r —t- 0: S ^ r^. The choice 
Ssing = —QSq IS the ouly one that satisfies Eq. f l28p in addition to the conditions (i) and (ii) 
above, for 7 = — |. 

While we believe that the above argument for the uniqueness of Ssmg = —GSq is quite 
plausible on physical grounds, the additional assumptions go beyond basic thermodynamic 
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Fi function of vq for different values of uq. Curve a corresponds 

to no = 0.01, curve b to = 0.1, curve c to no = 0.25, curve d to mq = up ~ 0.3861, curve e to 
uo = 0.44 and curve f to uq = uq ~ 0.4923. The maxima of curves d and f identify the points P 
and Q of Fig. 1, respectively. For up < uq < uq the curves have two, or more, local maxima, which 
correspond to different regular solutions with the same value of no. The maximum with the larger 
value of Vq is the global maximum. 




Figure 6: Stotiuo,vo, R)/R^^'^ as a function of logig vq for different values of uq. Curve A corresponds 
to Uq = 0.1, curve B to uq = 0.3, curve C to uq = 0.44 and curve D to uq = 0.75. 
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principles, and they may not provide a unique result when applied to other self-gravitating 
systems. It is necessary to perform a similar analysis on other systems before we can make a 
definite statement on this issue. 

In any case, the choice Ssing = —QSq is the simplest possible expression for the singularity 
entropy. The existence of a simple solution to a strongly constrained problem is an indication 
that the singularity entropy is not an artificial construct and that it represents a genuine 
feature of gravitational thermodynamics, that is fundamentally embedded into the structure 
of Einstein's equations. 

3.2 The physical interpretation of singularity entropy 

It is important to emphasize that the assignment of entropy to the singularities at r = 
does not necessarily imply an interpretation of these singularities as actual spacetime points. 
Such an assumption does not enter into the arguments employed for the identification of Ssing- 
In particular, our results do not require a commitment to the physical existence of naked 
singularities (however benign those may be), or to discount the possibility that quantum 
effects might smoothen out the singularity at the Planck scale. The key point here is that 
the term Sq in Eq. ( 125|) (and consequently Ssing) is associated to an interior boundary of the 
singular solutions to the Oppenheimer-Volkoff equations, and the application of the maximum 
entropy principle is well-defined whether the singularity corresponds to a spacetime point, or 
not. 

In order to explain this point, we consider that the interior solutions to the Oppenheimer- 
Volkoff equations are defined for values of r G [e, R], where e is a length cut-off. For example, 
e may correspond to a scale where quantum or statistical effects beyond relativity become 
important. Let us denote the corresponding spacetime manifold as M^. This can be embedded 
isometrically into a larger spacetime manifold M, which extends into the region r < e. For 
example, M may be the union of the spherically symmetric solution with a homogeneous 
solution with suitable junction conditions on the surface r = e — similar to a construction 
in Ref. [6j. The continued solution needs not contain any spacetime singularities. It is 
presumably determined by a microscopic theory of gravity, and its precise form is irrelevant 
to the arguments in this paper. 

The reason is that the maximum entropy principle applies, irrespective of the nature of 
the spacetime extension. The radiation entropy Srad is a sum of the term 5*1 defined on the 
exterior boundary r = R and of a term S^ defined on the interior boundary r = e. Since e 
can be taken arbitrarily close to zero, Se is well approximated by 5*0. In order to apply the 
maximum entropy principle, we consider a total entropy Stot = Srad + Sbound, where the term 
Sbound is defined on the boundary r = e and corresponds to the total entropy associated to the 
region r < e. The entropy Sbound is specified by the requirement that the maximum entropy 
principle holds. That is, we seek an expression for Sbound, such that Stot is maximized for fixed 
M and R. This is a well-defined maximization problem. Its solution is independent of the 
form of the extension in the region r < e. Whatever the physical configuration in the region 
r < e might be, its entropy equals the value Sbound that maximizes Stot- 

The solution to the maximization problem depends only on the functional form of Srad- 
Taking e arbitrarily close to 0, continuity implies that Srad is well approximated by the ex- 
pression fl23l) . which depends on the fields' values at the singularity r = 0. It follows that 
at the limit e — )■ 0, Sbound coincides with Ssing = — 65*0, as shown in Sec. 3.1. Thus, the 
limiting behavior of the fields near the singularity r = specifies the entropy Sging, even if 
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the singularity is not viewed as part of the spacetime. At present, the physical origin of Ssmg 
can only be a matter of conjecture: we can only identify it as a contribution to the entropy, 
defined on internal boundaries of solutions to Einstein' s equation, that is necessary for the 
existence of a consistent thermodynamic description. 

3.3 The implementation of the maximum- entropy principle 

The maximum-entropy principle asserts that equilibrium configurations correspond to the 
maximum value of entropy for a given value of the internal energy. If we identify the internal 
energy with the ADM mass M, namely, the energy measured by observers at asymptotic 
infinity, then the maximum-entropy principle applied to the entropy function Stot Eq- ( l27|) 
clearly implies that the regular solutions are the equilibrium configurations. 

The ADM mass M is a directly observable quantity for an observer outside the box. 
However, its identification with the thermodynamic internal energy is not mandatory. Even 
in non-gravitating systems, the internal energy does not always coincide with the total energy. 
Total energy depends on the reference frame: for equilibrium systems one usually identifies 
the internal energy with the total energy at the rest frame of the system's center of mass. In 
non-equilibrium systems, the internal energy is obtained from the total energy by subtracting 
the contribution of energy related to the dynamics of the flow |29j . 

A key problem in the identification of the ADM mass M with the internal energy is that 
the singular spacetimes considered here possess an internal boundary at the singularity r = 0. 
The mass function m(r) starts from a negative value — Mq at r = 0. Hence, one would expect 
that the value m{R) = M at the bounding box underestimates the internal energy of the 
system. 

Since we deal with static solutions to Einstein's equations, a better candidate for the 
internal energy U is the Komar integral, 

U = -^ [ dS,,^d^C, (29) 

where the integration is over all boundaries of a Cauchy surface E, dSf^i, is the area two- form 
on the boundary and ^ : = ^ . A general spherically symmetric metric is of the form 

ds' = -f{r)de + + r\de^ + sin^ Qdcj^^), (30) 

r 

where positive- valued function /(r) and m(r) the mass function. Evaluating the Komar 
integral for this class of metrics, we obtain 



U 



2 V r ^ 



(31) 



At r — )■ OO (outside the box), the metric is a Schwarzschild solution with ADM mass M. For r 
close to 0, the metric is given by Eq.([T9D. Hence, /(r) = (1 - Mo)/f f as r ^ 0. Evaluating 
Eq. (13T]) . we obtain 



U = M + ]^R ' vT^fY^^ = M + ^^oToo, (32) 
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where Too = Ty^l — uq is the temperature as measured by an observer at the asymptotic 
infinity and 5*0 is given by Eq. (l25ll . 

The maximum-entropy principle then imphes that the equihbrium configurations are de- 
termined by the maximization of the total entropy Stot for fixed values of the radius R and 
the Komar mass U. The scaling symmetry of Stot implies that the quantity 

a{uo, vo) = Stot{uo, vo, R)/R^^^, (33) 

does not depend on R. It is convenient to define the variable 

w := 2U/R = Uo + ^ (^) (ToVl^o, (34) 

where ctq = Ssing/R^^"^ depends only on uq and vq. 

Then, the equilibrium configurations correspond to pairs {uq^vq) that maximize (T(Mo,fo), 
subject to the constraint that w is constant. We denote the entropy-maximizing values of uq 
and Vq at fixed w as Ueqiw) and Veqiw) respectively. The entropy of the equihbrium states is 

S,,{U,R) = R^/^a[ueg{2U/R),v,g{2U/R)]. (35) 

For values of Mq such that regular solutions exist {uq G [0,uq]), entropy maximization at 

constant Wq gives the same result as entropy maximization at constant Uq. The proof of this 

statement is the following. 

We define the function amax{uo) '■= sup^^^ a{uo,vo), corresponding to the maximum of cr at 

fixed values of uq. As shown in Sec. 3.1, amax is an increasing function of uo, for uq G [0,mq]. 

Then, amax{uo) > (Tmax(Mo) for any u'q < Mo- 
Next, we consider a curve of constant w. Eq. (IMl) implies that for any point {uo,Vo) in 

this curve, Uq < w. Hence, for all points (mo^'^o) in the curve of constant w 

(^max (^ 

). (36) 

For w G [0,Mq], the maxima at constant Uq correspond to regular solutions, and equality in 
Eq. dSS]) is achieved for the regular solution that maximizes a at fixed value of Uq = w. Thus, 
for w G [0,mq], the maxima of a at constant uq and the maxima of a at constant w coincide. 

We conclude that entropy maximization at constant internal energy U and radius R iden- 
tifies the regular solutions as the equilibrium configurations. We note that the proof above 
does not require the explicit form of U, only the fact that U > M. 

One might inquire whether the identification of the Komar mass with the internal energy 
U removes the need for the introduction of the term Ssmg for the entropy of the singularities. 
Namely, whether the radiation entropy Srad has a maximum at the regular solutions for fixed 
values of U. The answer is negative: regular solutions correspond to local maxima of Srad at 
fixed U, but not to global maxima. Numerical evaluation shows that Srad — > cxd as fo — )■ 0, for 
any fixed value of U. The correct global behavior of the total entropy can only be guaranteed 
by the addition of the term Sgmg for the entropy of the singularities. 

To summarize, we argued that the internal energy of the singular solution should be iden- 
tified with the Komar mass Eq. fl5^ rather than the ADM mass M, and that the maxima of 
the total entropy Stot at fixed values of internal energy U and radius R correspond to regular 
solutions. Next, we proceed to identify the entropy-maximizing configurations in the sector 
of the thermodynamic state space where regular solutions do not exist. 
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Figure 7: The curves t^o = Veq{w) (left) and uq = Ueq{w) (right) corresponding to the equihbrium 
solutions. 

3.4 Equilibrium singular solutions 

In Fig. 7, we plot the functions Veqiw) and Ueq{w) that correspond to the maxima of a at 
fixed values of w. These functions determine the temperature T at the bounding box and the 
ADM mass of the equilibrium configurations for fixed internal energy U and box area AttR"^, 
as follows. 



We note that the equilibrium solutions span the whole range of parameters U and R. Moreover, 
there exists a unique solution for each pair (M, R) and at least one solution for each pair (T, R). 

As shown in Fig. 7, there exist three different 'phases' of equilibrium solutions. 

The first phase (R) corresponds to regular solutions for which Ueqiw) = w. It coincides 
with to the segment from the origin to the point Q in the curve of Fig.l. 

The second phase (S) corresponds to singular solutions with uq € {uq, 1) which maximize 
the entropy Stot for w G [uq, 1.2). It is interesting to note that these solutions also maximize 
Stot for fixed uq. Hence, the same set of equilibrium solutions are obtained even by identifying 
the internal energy with the ADM mass M. The S-phase consists of type II solutions , 
except for values of Uq near 1, where the qualitative distinction between types I and II is not 
meaningful. 

The third phase (D) consists of degenerate configurations. These configurations are degen- 
erate in the sense that they all correspond to the same value of uq 1_ and are distinguished 
by different values of w > 1.2. Strictly speaking, a configuration of type D is not a solution of 
the Oppenheimer-Volkoff equations, because no static solution exists if Uq = 1. In a solution 
with Uq = 1, the bounding box would lie in a null (rather than a spacelike) hypersurface, 
thus violating the set-up of the system. The D-phase is a boundary of the thermodynamical 
state space that is obtained from solutions to the Oppenheimer-Volkoff equations through a 
limiting procedure. In practice, the introduction of a cut-off length-scale /, restricting the 
possible values of Uq to be smaller than 1 — 21/R, allows for an identification of the global 
maxima of entropy with the solutions corresponding to Mq = 1 — 21/R, at least for values of 
temperature T such that IT « 1. 



M{U,R) = -Rueq{2U/R). 




(37) 
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Figure 8: a^q is plotted as a function of w. The first derivative of a^q is discontinuous at the value 
ol w = uq corresponding to the transition from regular to singular solutions, but aeq remains an 
increasing function until w ~ 0.51. Then a^q decays to zero with increasing w. Degenerate solutions 
all have zero entropy. 

We believe that the degenerate solutions should be physically interpreted as a phase of 
complete gravitational collapse. The idealization of a box bounding self-gravitating radiation 
becomes unphysical, because for solutions at — ?■ 1 the stresses on the box tend to infinity. 
We do not obtain an actual black hole solution, because the equation of state of radiation does 
not admit solutions with vanishing density at finite radius, hence, it cannot describe compact 
objects, unbounded by a box. We expect that systems with different equations of state have 
a richer and more physically relevant phase structure at the corresponding limit. 

In Fig. 8, we plot the function (Jeq{w) = Seq/B?^"^-, where Seq{U, R) is the entropy of 
equilibrium states, defined by Eq. (135|) . For regular solutions, aeq is an increasing function 
of w; for singular solutions a^q remains an increasing function of w up to w ~ 0.51 where it 
starts decreasing and vanishes at w = 1.2. Degenerate solutions have vanishing entropy. 

Fig. 8 shows that for any given 'radius' R there is a maximal value of Stot that approxi- 
mately equals 0.823i?^/^. This contrasts the behavior of entropy in non-gravitating systems, 
where the entropy remains an increasing function of the internal energy at constant volume 
and has no maximum. On the other hand, the existence of a maximal value for the equilib- 
rium entropy at constant R is in accordance with the idea of an entropy bound for localized 
gravitating systems [30] . 

There is also a global minimum to Stot for fixed R, approximately equal to — 1.65-R^/^. Since 
we chose the scale of entropy such that Stot = for the solutions with M = 0, the positivity 
of entropy suggests a bound on the gravitational entropy Sfiat{R) of flat space within a box 
of radius R, namely, Sfiat{R) > 1.65i?'^/^. This bound is not particularly restrictive. For 
example, any expression for SfiatiR) of the form cR"', for 7 > 3/2 and c of order unity satisfies 
this bound for all radii R larger than the Planck scale. 

The decrease of entropy with increasing internal energy U at constant radius R, character- 
izing the S-phase, would be quite unphysical in ordinary thermodynamics, because it would 
imply a negative temperature. However, in a self-gravitating system the equilibrium configu- 
ration is inhomogeneous and there is no guarantee that the variable T, corresponding to the 
temperature at the bounding box, has a simple relation to the entropy. 

In the regular phase, there exists only a single boundary, and the relation between tem- 
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perature and entropy is formally analogous to the one in extensive thermodynamics, namely, 
T^^ = {dSeq/dU)R [6J. This simple relation is lost in the singular phase, because singular 
solutions are characterized by an internal boundary at the singularity. In fact, the singular 
solutions exhibit a feature very much unique to gravitational systems: the boundary of the 
system consists of two disconnected components, namely, the box at r = R and the singular- 
ity at r = 0. The two components of the boundary have different temperature, and yet the 
system is in equilibrium with no energy flow. For this reason, the relation between the tem- 
perature(s) of the system and the derivatives of entropy is expected to be much more complex. 
We have not been able to obtain a physically intuitive form of this relation. The study of 
other gravitating system with an inner boundary is necessary, in order to better understand 
this issue. 

The results obtained above for the non-equilibrium singular solutions rely on the identi- 
fication of the internal energy U with the Komar mass of the solutions. The Komar mass for 
radiation equals 2 y/—gd^xp, where S is a Cauchy surface, i.e., it includes only the contri- 
bution of radiation. In principle, it is possible that the singularity contributes an additional 
term Using defined from the values of the field variables at r = 0, analogous to the term Ssing 
in total entropy. 

We have found no argument that allows us to uniquely determine such a term. The reason 
is that the key thermodynamic properties are not significantly affected by the choice of U. 
In particular, the proof in Sec. 3.3, that the regular solutions maximize the total entropy for 
fixed U and R, holds for any choice of U > M. 

The simplest choice for Using is a term proportional to RTooSo, as in Eq. (1321) . This means 
that the internal energy is of the form 



for some free parameter A. We have verified numerically that the curves of Fig. 7 remain 
unchanged for a large range of values of A (roughly, for A e [0, 10]), i.e., the entropy maxima 
at constant U are insensitive to the choice of A, even for the singular solutions. Hence, the only 
effect of a different definition of f/ is a reparameterization of the variables w in the plots of Fig. 
7; in particular, the decrease of the equilibrium entropy in the S-phase remains unchanged. 

3.5 Phase transitions 

Next, we examine the transitions between the different phases R, S, and D of self-gravitating 
radiation. In systems with long-range forces, the usual arguments about the equivalence of 
the microcanonical and the canonical distribution (at the thermodynamic limit) fail to apply 
|16j . Hence, the phase structure of the equilibrium states of an isolated system may well differ 
from that of a system in contact with a thermal reservoir at constant temperature T. 

We first consider an isolated box containing self-gravitating radiation. The relevant ther- 
modynamic potential is the equilibrium entropy Seq- The phase structure is described by Fig. 
8. There, we note two transition points, one at w = uq corresponding to the transition from 
the R-phase to the S-phase, and another at w = 1.2 corresponding to the transition from the 
S- phase to the D-phase. Both phase transitions are continuous. 



U = M + -RT^So, 



(39) 
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Figure 9: The scaled Helmholtz free-energy / = F/E?/'^ is plotted as a function of the scaled local 
temperature TR^^^ for all equilibrium solutions. 

The first derivative of the temperature T a.t w = uq has an infinite discontinuity ^22j. We 
find that (|^)_ ~ -0.57R-^/^{uq - w)'^ as w ^ uq and (|^)^ ~ -0.28R-'^/^. Thus, there 
is a finite discontinuity of the heat capacity C = at w = uq: C+ — C- ^ —5.05/?'^/^. 
Numerical evaluation has not shown any discontinuities in the first derivatives of the basic 
thermodynamic variables at the second critical point w = 1.2. 

Since the singularity entropy Sging vanishes for the regular solutions and it is non-zero for 
the singular and degenerate solutions, it is natural to associate |S'si„g| with the order parameter 
for the system. We find that near the critical point w = uq, \Ssing\ ~ 1.73i?'^/^(w ~ uq). The 
behavior of the system near the critical point suggests a phase transition described by the 
critical exponents of mean-field theory (with \/\Ssing\ as the order parameter), as is to be 
expected from a treatment that relies on the maximum-entropy principle. 

Next, we consider the box in contact with a thermal reservoir at temperature T. As seen 
from Fig. 7, there exists a range of values of Vq, and hence a range of values of T, where 
all three phases coexist. The equilibrium phase corresponds to the global minimum of the 
Helmholtz free energy F = U — TooSeq- As in ordinary equilibrium thermodynamics, this 
follows from the entropy maximization of the total system including the box with radiation 
and the thermal reservoir, and taking into account that an exchange of energy SQ at the 
boundary corresponds to a change 6U = a/1 — Uo6Q of the internal energy of the box. In Fig. 
9, the scaled Helmholtz free energy / = F/R^/"^ is plotted as a function of the scaled local 
temperature 

The regular phase corresponds to a global minimum of the free energy for all values TR^^^ G 
[0,0,406). For higher temperatures only the degenerate phase exists. The critical point 
TR^/"^ ~ 0.406 corresponds to the point P of Fig. 1. The segment PQ of Fig. 1 as well 
as the S-phase do not correspond to global minima in the free energy, and thus they do not 
define stable equilibrium states. Hence, when the box is in contact with a thermal reservoir, 
there is a single critical point where the system undergoes a first- order phase transition, from 
the regular to the degenerate phase. The latent heat for this transition (as measured by an 
observer at infinity) is defined as L := Ud — Ur, where Ud and Ur are the internal energies 
at the critical point for the D-phase and the R-phase respectively. We find that L ~ 0.237i?. 
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The equilibrium states determined in Fig. 9 are characterized by non-negative values of 
^ 0). This is in marked contrast to the case of the isolated box considered earlier. 
This property reflects the positivity of the heat capacity characterizing systems described by 
the canonical distribution |8]. 

The phase structure of self-gravitating radiation is qualitatively similar to the phase struc- 
ture of isothermal Newtonian spheres, regularized at small scales — see, Ref. [10] for a review 
of such transitions using the microcanonical and canonical ensembles. We note that the ap- 
proach we employ in this paper provides only a rough approximation to the physics of the 
phase transition. The reason is that the maximum-entropy principle cannot account for the 
effect of thermodynamic fluctuations that are essential for the precise description of critical 
phenomena [20]. A more reliable approach for the study of phase transitions should involve 
statistical thermodynamics. The entropy functional Stot allows the definition of equilibrium 
states as averages with respect to the probability distribution exp[5'tof(f/, T, i?)]. Thus, the 
contribution of the fluctuations near the critical points would be incorporated into the ther- 
modynamic description. 

3.6 Relation to Bekenstein-Hawking entropy 

Black hole entropy has provided the most important indication of a fundamental relation 
between gravity and thermodynamics. Therefore, it is natural to inquire whether the entropy 
of singularities, proposed here, is related to the Bekenstein-Hawking entropy. 

The arguments developed in this paper involve only the basic principles of general relativity 
and thermodynamics. At no point has it been necessary to introduce quantum theory, and, 
consequently h does not appear in any of the expressions. (It only appears implicitly through 
the Stefan-Boltzmann constant.) Hence, the Bekenstein-Hawking entropy cannot appear in 
the present context, because it depends explicitly on The singularities of self-gravitating 
radiation have less entropy than locally Minkowski an spacetime, while black hole horizons 
have much larger entropy. Presumably, they originate from different sectors in the space of 
gravitational microstates. 

In Sec. 2.3.1, we noted that the system admits a class of trivial solutions for uq G (0, 1) 
and f = that correspond to a Schwarzschild horizon of radius rs = 2M inside the box and 
no radiation. These solutions form another, fourth, phase of the system. Within the classical 
theory, their entropy is zero; hence, by the maximum-entropy principle they do not define 
equilibrium solutions. However, one might inquire whether an assignment of Bekenstein- 
Hawking entropy Sbh = 7rr|/^ to these solutions can be consistently incorporated into the 
present framework. The answer is negative for the following reasons. 

A black hole horizon, endowed with thermodynamic properties, is expected to emit Hawk- 
ing radiation. The presence of a bounding box implies that this radiation does not escape to 
infinity. Thus, the equilibrium configuration should correspond to the horizon coexisting with 
its Hawking radiation. The local temperature T at the bounding box would equal the Hawking 
temperature blueshifted by a factor of \/l — Uq, as in Ref. [31]. But, as we showed in Sec. 2.3, 
there is no solution of the Oppenheimer-Volkoff equations with both non-zero temperature 
at the boundary and an horizon. We would have to modify the constitutive equations of the 
system in order to introduce the Bekenstein-Hawking entropy into the present set-up. Such 
modifications could arise from the consideration of quantum corrections to Einstein's equa- 
tions, or from the inclusion of the effect of acceleration radiation, or from modifications to the 
equation of state for radiation near the horizon [32]. However, any modification would affect 
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the entropy of the system (both of radiation and of the singularities) , possibly introducing ad- 
ditional terms with different scaling properties. As a matter of fact, such considerations might 
provide an important generalization of the results presented here, namely, the determination 
of constitutive equations that incorporate the Bekenstein- Hawking entropy into the classical 
thermodynamic description. 

4 Conclusions 

We have argued that the requirement of thermodynamical consistency, when applied to self- 
gravitating systems, can provide novel information about the thermodynamic properties of 
the gravitational field. We applied this idea to a concrete system, namely, static, spherically 
symmetric solutions to Einstein's equations describing self-gravitating radiation in a box. 
We showed that the only way to describe this system in a thermodynamically consistent way 
requires the assignment a specific expression Ssing to the entropy of the spacetime singularities. 
While this result is restricted to a specific class of self-gravitating system, it does demonstrate 
a point of fundamental importance: that horizons are not the only entities of gravity theory 
to which entropy can be assigned. This implies the intrinsic connection between gravity and 
thermodynamics is not restricted in the black hole context. 

Of equal importance is the fact that the expression Ssing for the singularity entropy is 
identified through stringent constraints of physical and mathematical consistency. The exis- 
tence of a functional satisfying such constraints is a highly non-trivial result. It suggests that 
the thermodynamic description of the gravitational field is fundamentally embedded in Ein- 
stein's equations, and the rationale developed here can be applied to other physical systems. 
In particular, we expect that the justification of the regularity condition through maximum- 
entropy principle holds in general spherically symmetric solutions to Einstein's equations, as 
long as the equations of state are thermodynamically consistent. The study of a large class 
of such systems will enable the solidification of this hypothesis. The eventual aim is to iden- 
tify gravitational entropy within the broadest framework where the principles of equilibrium 
thermodynamics are applicable, namely, for general static solutions to Einstein's equations. 

The main limitation of the present results is the fact that they are based on numerical 
solution of the Oppenheimer-VolkofT equations, rather than a full analytic derivation. The 
application of techniques from the theory of differential equations is an important requirement 
of further research, because it can provide a firmer analytic control on the properties of singular 
solutions. The definition of the singularity entropy requires very specific behavior of the fields 
near the singularity. It is essential to understand the mathematical origins of this behavior, 
in order to proceed towards directions of increasing generality. 

It is important to emphasize that the methodology developed here employs only assump- 
tions from classical theories, namely, general relativity and thermodynamics. Thus, any results 
obtained from this method can, in principle, be employed in order to constrain candidate quan- 
tum theories of gravity. The study of phase transitions is perhaps the most important in this 
direction. As explained in Sec. 3.5, the maximum-entropy principle is not adequate for the 
description of the systems near the critical point. However, the definition of a consistent 
entropy functional provides the basis for a treatment of the critical behavior using statistical 
thermodynamics. If the gravitational phase transitions exhibit universality, then the identi- 
fication of a universality class could provide novel information about the symmetries of the 
underlying microscopic theory. 
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This paper is partly motivated by a change of perspective towards singularities. We do not 
view them as pathologies of classical general relativity, to be resolved by a quantum theory of 
gravity, but as indications of a different thermodynamic phase for the gravitational field, which 
appears when the external constraints are not compatible with locally Minkowskian geometry. 
This perspective might have implications in the theory of stellar stability, where physical 
solutions are standardly selected by the regularity assumption. A thermodynamic treatment 
of systems with realistic equations of state, could uncover a new phase of singular equilibrium 
solutions, similar to the S-phase appearing in self-gravitating radiation. If such a phase turns 
out to be thermodynamically stable when the system is in contact with heat and/or pressure 
reservoirs, it would suggest the existence of a new type of compact astrophysical objects, that 
interpolates between ordinary stars (regular solutions to Einstein's equation) and black holes. 
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